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Surface plasmon-polariton (SPP) excitations of metal-dielectric interfaces are a fundamental light-matter 
interaction which has attracted interest as a route to spatial confinement of light far beyond that offered by 
conventional dielectric optical devices. Conventionally, SPPs have been studied in noble-metal structures, 
where the SPPs are intrinsically bound to a 2D metal-dielectric interface. Meanwhile, recent advances in the 
growth of hybrid 2D crystals, which comprise laterally connected domains of distinct atomically thin 
materials, provide the first realistic platform on which a 2D metal-dielectric system with a truly ID 
metal-dielectric interface can be achieved. Here we show for the first time that ID metal-dielectric interfaces 
support a fundamental ID plasmonic mode (IDSPP) which exhibits cutoff behavior that provides 
dramatically improved light confinement in 2D systems. The IDSPP constitutes a new basic category of 
plasmon as the missing ID member of the plasmon family: 3D bulk plasmon, 2DSPP, IDSPP, and OD 
localized SP. 

Low- dimensional collective electron excitations at metal surfaces, the so-called surface plasmon-polaritons 
(SPPs), provide a route towards tremendous electric field enhancement^ and spatial confmement\ which can 
dramatically enhance light-matter interactions^. Conventionally, SPPs have been studied in low-dimensional 
noble-metal structures, where the SPP is intrinsically bound to a 2D metal-dielectric interface\ Although recently, 
graphene^, the atomically thin 2D hexagonal crystal of carbon atoms, has rapidly gained interest as a long-sought- 
after plasmonic material^'^^ alternative to noble metals due to initial predictions^ of exceptional electric field 
confinement of the intrinsic 2D plasmonic excitations of doped graphene sheets^"^. Quasi-low-dimensional 
schemes aiming to further increase the electric field confinement of graphene plasmons, such as excitations of 
graphene ribbon structures^"^^, p-n junctions^^'^^, discs^'^^'^^ and nanoresonators^^, have been suggested. 

Meanwhile, at the cutting edge of materials science are hybrid 2D crystals^^"^^ comprising laterally connected 
domains of distinct atomically thin materials, and which are mechanically continuous over macroscopic 
domains. Moreover, studies have shown that compositional transition can occur over atomic scale distances^^"^^ 
at a crystalline junction^^ connecting the neighbouring domains. Efforts in this direction have been largely driven 
by prospects to further complement the outstanding properties of graphene^^, and towards its implementation 
with other 2D materials such as the insulating hexagonal boron nitride (h-BN), for atomically thin electrical 
circuitry^^'^^. However, yet to be recognized is that hybrid 2D crystals further provide the first realistic platform on 
which a 2D metal-dielectric system with truly ID metal -dielectric (IDMD) interface - that is, across which the 
sign of imaginary part of conductivity changes from positive (i.e., metallic response), to negative (i.e., dielectric 
response) - could be achieved and probed at optical frequencies. Although, while such a ID interface is of 
remarkable fundamental simplicity, constituting the low-dimensional counterpart of the bulk metal-dielectric 
interface, the existence and nature of plasmonic excitations at IDMD interfaces remain unknown. 

In this study, we show that IDMD interfaces in general 2D metal-dielectric systems support a fundamental and 
unique ID plasmonic mode (IDSPP). Through an illustrative example on a hybrid graphene/graphene platform, 
we show how unique cutoff behaviour of IDSPPs could allow for a dramatic improvement in the electric field 
confinement of plasmons in 2D systems, exceeding that of previously predicted quasi- ID and ID plasmonic 
excitations in graphene^"^^. In the considered example, we demonstrate confinement of the electric field intensity 
to modal areas over one million times smaller than the diffraction limit. Furthermore, from a fundamental 
perspective, the IDSPP, as the fundamental excitation of a IDMD interface, constitutes a new basic category 
of plasmon on its own: the missing ID member of the plasmon family (3D bulk plasmon, 2DSPP, IDSPP, and OD 
localised SP). 
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Results 

We start by considering a 2D metal -dielectric (MD) system in they = 
0 plane consisting of two semi-infmite domains that are laterally 
connected along the z axis with sheet conductivities of a^^^ = g 
+ /a^^)" (x < 0) and a^^^ = a^^^' + ia^^^" (x > 0) (see inset in Fig. la) 
and immersed in a uniform dielectric with relative dielectric permit- 
tivity 8. Without loss of generality, we take regions L and R as the 
metallic (g^^^" > 0) and dielectric (g^^^" < 0) domains, respectively. At 
this point, to simplify our discussion, we omit ohmic and interband 
losses in both domains (a^^'^^' = 0; the lossy case is considered below 
and in the Supplementary Information). According to the theory of 
Volkov and Mikhailov^^'^^, the dispersion relation of plasmons prop- 
agating along the IDMD interface {z = 0) in the quasi-static limit can 
be shown to be (see Methods): 

deinl =/2m-l7rV2, 1 

where o"^^^" = Ka^^^" {K < 0), the normalised effective index N = 
nln2D {N > 0; n = q/ko; = w/c) is the effective index of the ID 
interface plasmon with wavenumber q normalised to that of the 
planar 2D transverse magnetic (TM) plasmons of region L, /I2D = 
IceeqIg^^^" (Ref 5), and m is an integer. 

While equation (1) lacks a closed-form solution, our numerical 
solution indicates (see Methods) a dispersion relation of the form 

q^2^3j2EEQOjl{G^^^" + G^^^"^ , and N{K) is plotted in Fig. la. We 

find that non-leaky plasmons (i.e., AT is purely real) exist strictly in the 
window 0< —K< l.AtiC = 0, the plasmon wavenumber takes its 

minimum value as that of the bare edge plasmon^^'^^ q^ ~ ^3/2^20^0 • 
As ^ 1, that is, when the electric susceptibilities on the two sides 
of the interface become equal (|o-^^^"| g^^^"), the plasmon wave 
number diverges. Interestingly, this divergent behaviour reveals a 
strong connection to the conventional 2D surface plasmon polari- 
tons (SPPs) locaUsed to bulk metal-dielectric interfaces'^. Indeed, in 
the absence of ohmic loss, the wavenumber of the SPPs (<^spp = 
^spp^o) diverges at the surface plasmon (SP) frequency Wsp given 
by the non-retarded SP condition'^: + £m(^^sp) = 0;8d> 0 and 8^ 
< 0 are the relative permittivities of the dielectric and the metal, 
respectively. Likewise, at the IDMD interface, we note the cutoff 
condition K = - 1, or g^^^" + a^^^" = 0 (a^^'^^' = 0). This reveals that 
the plasmonic mode of a IDMD interface is a ID manifestation of 
SPPs in 2D systems (IDSPP). An important consequence of the 
reduced dimensionality can be observed in the respective divergence 



behaviours: the SPP behaviour takes the form ngpp ^ 1 / a/As (As = 8^ 
~ |£m|)j which diverges much less rapidly than that of the IDSPP, 

n~l/A(7(A(7 = (7^^)" - \g^^^"\). 

We note that in the earlier works of Volkov and Mikhailov'^'^^ 
they considered only conventional 2D electron gases (2DEG) with 
G > 0 i.e., where both domains (L and R) are metallic. In that 
case, while there is no solution to equation (1) for purely real N, there 
is a solution for complex N which corresponds to a so-called inter- 
edge plasmon that is damped as it leaks into the extended 2D plas- 
mons supported by the domain with the smaller (positive) imaginary 
part of conductivity'^. In contrast, the IDSPP (existing only in the 2D 
metal-dielectric system when 0 < —K< 1) is entirely non-leaky, and 
completely bound to the IDMD interface. It is also important to note 
a clear distinction between the IDSPP, and the recently predicted 
quasi- ID plasmonic excitations of graphene p-n junctions^^'^^, 
according to fundamental differences between their respective phys- 
ical origins. The multiple plasmon modes of a graphene p-n junc- 
^jQj^i4,i5 their manifestation to a spatial gradient of the imaginary 
part of conductivity (i.e., g"), much like the multiple modes of a 
graded-index waveguide, and the graphene is everywhere metallic. 
In contrast, the IDSPP manifests exclusively as a consequence of the 
metallic-to-dielectric transition across the IDMD interface, either 
side of which the conductivity is spatially uniform. In this regard, 
we predict that the IR edge wave observed graphically in Vakil et al. 
(Ref 9) is in fact due to the excitation of IDSPPs, rather than p-n 
junction plasmons as suggested therein. 

The lossy IDSPP dispersion (i.e., when either g^^^', g^^^' 0) is 
readily obtained on introducing the complex parameter K defined in 
terms of the complex conductivities K = g^^^Ig^^^ = K' + iK", and 
solving for the complex normalised effective index N = nln2r> = N' 
+ iN" (where = i2c8o/G^^^ - Ref. 5). In the limit that the loss is 
wholly confined to the dielectric domain (g^^^' = 0, g^^^' 0), we 
obtain the convenient and meaningful expressions K' = a^^^" j q^^^" 
(i.e., same as the lossless case), and K' = —g^^^'Ig^^^". Thus we may 
imagine that fixing K' and simultaneously increasing \K'\ corre- 
sponds to fixing of the imaginary part of conductivity of the L and 
R domains, and linearly increasing the loss {g^^^') in the dielectric 
domain. Conversely, fixing \K'\ and increasing K' is equivalent to 
fixing G^^^' and linearly increasing cr^^^". Meanwhile, the quantity 
N'/(2nN") reduces to the propagation length defined as the 
number of optical cycles within one exponential decay length 
(l=Re(<5)/[27iIm(<^)]). We note that this physical scenario is particu- 
larly relevant to a hybrid graphene/graphene metal- dielectric system 
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Figure 1 | IDMD interface in a hybrid 2D system and dispersion of IDSPPs. (a), Effective index n of the IDMD interface plasmon ( IDSPP) normalised 
to that of planar 2D plasmons in the metallic domain L (/72d) as a function of = q^^^" j g^^^" (ct^l)" > Q; o^^^" < 0). Inset: Schematic of the considered 2D 
metal-dielectric system comprising two semi-infinite domains L (x < 0; metallic) and R (x > 0; dielectric) laterally connected along the z axis, with 
sheet conductivities of o^^^ and o^^\ respectively. We consider a plasmon with a wavevector q = that propagates along and is localised to the ID 
interface (coinciding with the zaxis). (b), K' dependence of N'/{2nN") at indicated values of K'. Inset: lossy dispersion of N' for example values of loss K" 
= —0.06 (green), K" = —0.12 (red), as compared to lossless dispersion (blue). 



SCIENTIFIC REPORTS | 4 : 4536 | DOI: 1 0.1 038/srep04536 



2 



at finite temperature (see Supplementary Information) where the 
dielectric response of domain R is accompanied by interband losses 
which strongly dominate the overall loss (i.e., relative to ohmic loss in 
either domain). 

The green and red curves in the inset of Fig. lb show N' when loss 
is introduced through two example values K" = —0.06 and K" = 
— 0.12. Similarly to the behavior of conventional SPPs on a lossy 
metal surface^^, the loss is shown to destroy the divergent cutoff 
behavior of the IDSPP owing to heavy damping, the onset of which 
occurs in some vicinity of K' that is specific to K" (compare where N' 
begins to decrease at the two values of K"). To illustrate and quantify 
the onset of IDSPP damping, the main window of Fig. lb shows the 
K" dependence of the IDSPP propagation length at several fixed 
values of K'.Ata given value of K", the propagation length decreases 
with an increase in K\ which is typical of plasmons as their wave 
number increases (recalling that q increases with K'). We note that, at 
each value of K\ the curve is cutoff at the x-axis corresponding to the 
propagation length of just one wavelength (i.e., a strongly damped 
IDSPP). It is illustrative to compare the K" cutoff with a particular K' 

e.g., when K' = —0.6, the cutoff is K" 0.06, corresponding to K"/ 

K' = —g^^^'/g^^^" — 0.1. Obviously, the loss is more forgiving for 
smaller K\ e.g. when K' = —0.2, the cutoff is X" 0.12, corres- 
ponding to -a^^^'/a^^^" ~ 0.6. 

Now we present an illustrative and topical example of IDSPPs in a 
hybrid graphene/graphene metal- dielectric system, in which the ran- 
dom phase approximation (RPA) provides convenient closed-form 
expressions for the optical conductivity in both domains^"^'^^. The 
IDMD interface is achieved by appropriate doping of both graphene 
domains; in the local and zero-temperature limit considered 
here (see Ref 30 and Methods), the doping would correspond to 
Id^^Vhw > 0.6 and 0.5 < ^^^V^zco < 0.6, where ^^^'^^ is the Fermi 
energy of domain L or R, respectively, and details at finite tempera- 
tures are given in the Supplementary Information. Our numerical 
calculations suggest that the existence of IDSPPs is not critically 
dependent on an abrupt sign change of g" at the interface (see 
Supplementary Information), so hybrid graphene systems that sup- 
port IDSPPs may also be realised by proposed electrical gating^ or 
substrate controlled^^ schemes for spatially non-uniform graphene 
doping in addition to patterned growth^^"^^. We note that in an 
electrical gating scheme^, the absence of a well-defined edge separat- 
ing the metallic and dielectric graphene domains would reduce edge 
effects such as those associated with zigzag edges in graphene and 
known to be responsible for strong plasmon damping (see for 
example Ref 32). 

The conductivities on either side of the interface are completely 
determined by the respective normalised frequencies^" Q^^'^^ = hw/ 
^(L,R) j^^g been shown that the dispersion of IDSPPs (relative to 
^2d) is completely determined by the dimensionless conductivity K, 
this behaviour suggests a flexible tunability of n through the variation 
of either a^^^" (fi^^^) or a^^^" (fi'-^^) with the doping of either graphene 
domain. We demonstrate this tunability and simultaneously verify 
our quasi-static analytical results by employing a mode-solver tool in 
the numerical finite element method (FEM) package COMSOL, 
which includes retardation effects. In Fig. 2a, we show the calculated 
dependence of n on /n^^^ for several fixed values of ^^^^ (see respective 
curves) and the spatial evolution of the IDSPP electric field compo- 
nents Ex and Ey. Note the normalisation of fi^^'^^ in terms of Hw; this 
was verified by comparing all results at the frequencies / = 20 THz 
and 80 THz. One can observe excellent agreement between the dis- 
persion analytically predicted by equation (1) (solid curves) and the 
numerically determined data points. The monotonic increase in n as 
fi^^^ decreases (at fixed ^^^0 can be readily understood in terms of the 
corresponding decrease in a^^^^and thus Act (recalling that n ~ 1/Aa). 
Each curve diverges asymptotically towards a cutoff value of /n^^^ 
given by the condition a^^\fi^^^) + a^^\fi^^^) = 0. The cutoff value 
of fi^^^ approaches 0.5hw as fi^^^ increases. When the plasmon energy 
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Figure 2 | IDSPPs in a graphene/graphene hybrid 2D system, (a), 
Effective index n = q/ko of the IDSPP on a sample hybrid graphene/ 
graphene system as a function of (Q^^^)"^ = ji^^yhco; the different curves 
correspond to (Q^^^)"^ = ii^^^/Hw = 2, 3, 4, and 5. Excellent agreement is 
observed between the analytic dispersion equation (1) (solid curves) and 
FEM simulations (data points). Left inset: Spatial evolution of the IDSPP 
electric field components (xz slice) and Ey {yz slice); (Q^^^)~^ = 2 and 
(Q(^))~i ~ 0.5. Right inset: Dependence of the net power flow P (along z) 
on the effective index n of the IDSPP, normalised to that of the graphene 
bare edge plasmon (Pe> ^e)- (b), Electric field components of the IDSPP for 
the indicated parameters; the profiles on the left correspond to the IDSPP 
near the cutoff, while those on the right correspond to the graphene bare 
edge plasmon (GEP) [the horizontal line indicates the graphene plane; the 
red (blue) colour denotes a positive (negative) sign] . 

hoj and the Fermi energy p!^^ are of approximately the same order of 
magnitude (corresponding to the typical practical situation), the 
cutoff is very close to 0.5 because of the relative strengths of the 
interband and intraband terms in the graphene conductivity (see 
Methods). 

The electric field profiles in Fig. 2b combined with the plot of the 
net power flow (inset in Fig. 2a; normalised to the graphene bare edge 
plasmon^^ (GEP), i.e., when g^^^ = 0) demonstrate the nature of the 
cutoff dynamics at the IDMD interface. The net power flow parallel 
to the interface is given by the integral of Sz = E^Hy^ — EyH^^ over 
the xy plane. is symmetric about the x axis, although the terms 
E^Hy"^ and EyH^"^ both exhibit odd sign parity about the y axis; thus, 
cumulative integration on either side of the interface leads to their 
respective partial cancellation. This cancellation is weak in the case of 
GEPs because of the strong field asymmetry, and therefore, the net 
power flow is significant. Gonversely, as the cutoff condition {g^^^" 
— G^^^") is approached, the symmetrisation of the field components 
caused by equal charge screening on either side of the interface 
(analogous to conventional SPPs near the surface plasmon fre- 
quency^^) leads to a strong reduction in the net power flow as the 
group velocity approaches zero. 
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Figure 3 shows the evolution of the IDSPP mode cross section (see 
caption for details) as the cutoff is approached, in comparison with 
the GEP. We note the tremendous locaUsation of the electric field 
intensity (lower inset) to regions from ~ 5 X IO'^Xq^ to A ~ 6 X 
lO~'^Xo^ as \K\ increases in the range 0 < < 0.87. Thus, the spatial 
confinement of the electric field intensity of the IDSPP in the hybrid 
graphene/graphene system is on the order of one million times smal- 
ler than the diffraction limit (~^o^; is the vacuum wavelength); a 
two-order-of-magnitude improvement over the graphene bare edge 
plasmon (i.e., when K = 0). 

Discussion 

We emphasise that the above results are not limited to the presented 
example of graphene but are generalizable to any isolated IDMD 
interface in a hybrid 2D system characterised by the same value of 
K. IDMD interfaces in graphene/graphene hybrid 2D systems at 
finite temperature are intrinsically lossy because of the presence of 
interband transitions, which are essential to achieving the dielectric 
character in one domain (i.e., cr^^^" < 0; iC < 0), and further exhibit a 
temperature -dependent upper limit on (and \K\). Neverthe- 

less, we predict the propagation lengths of IDSPPs on hybrid 
graphene systems to be several plasmon wavelengths at room tem- 
perature (T = 300 K) and —100 plasmon wavelengths at the tem- 
perature of Hquid nitrogen {T = 80 K); for example, at a free space 
wavelength X = 1.8 |im we find / -lOOnm {fi^^^ = 1 eV, fi^^^ = 
0.42 eV, T = 300 K) and / ~ 2.6 |im (^^^^ = 1 eV, fi^^^ = 0.40 eV, 
T = 80 K), respectively (see Supplementary Information for full 
details). On the other hand, graphene/h-BN hybrid structures are 
promising systems for supporting IDSPPs with large propagation 
lengths. Because of the dielectric character of h-BN, one could realise 
a graphene/h-BN IDMD interface with separate tunability of the 
graphene conductivity in real time using external gating, and the 
plasmon propagation length in this system would be limited only 




Figure 3 | Enhanced modal confinement of IDSPPs. Mode cross section 
of IDSPPs as the cutoff condition is approached (blue-dashed K= 0^ red 
K = —0.87). Concentric contours, from the outermost to the innermost 
contour, correspond to (Q^^^"' = 0.6, 0.58, 0.56, 0.54, 0.52, 0.501, and 
0.5002. The blue-dashed contour exactly coincides with that of the 
graphene bare edge plasmon (GEP). Each contour is defined as a line 
lEP = constant enclosing an area A (in the xy plane), such that 



||E|^dS/| 



|E|^<iS = 0.8. Insets: (top) Mode cross section area A 



normalised to that of the GEP Aq as a function of (Q^^^) ^ and (bottom) 
normalised intensity of the electric field (2 nm above the graphene) for the 
extreme cases (Q^^^)"' = 0.6 (blue-dashed) and (Q^^^"^ = 0.5002 (red). 
/= 50 THz; (Q(L))-^ = 2. 



by the small ohmic losses in high-mobility graphene. Although in the 
case of an atomically sharp transition between neighbouring 
domains in hybrid 2D crystals, consideration of edge related effects 
(such as electronic edge states^^) would be important to accurately 
determine the IDSPP dispersion and loss properties. 

In summary, we have predicted the existence of a fundamental ID 
plasmonic mode of ID metal-dielectric interfaces in 2D systems 
(IDSPP). The effective index of IDSPPs diverges asymptotically 
towards a cutoff as the magnitudes of the susceptibilities on the 
two sides of the interface become equal, in striking analogy to 
conventional SPPs at bulk metal-dielectric interfaces. On a sample 
2D metal-dielectric graphene/graphene platform, highly sensitive 
tunability of the IDSPP dispersion was demonstrated via the doping 
of either graphene domain, which further allowed for the achieve- 
ment of spatial confinement of electric-field intensity to regions 
orders of magnitude smaller than that of the plasmonic excitations 
of a bare graphene edge (GEP). The unique and tuneable cutoff 
behaviour of IDSPPs thus presents a means for dramatic enhance- 
ment of light confinement in 2D systems over other ID or quasi- ID 
excitations, such as GEPs or plasmonic excitations of graphene p-n 
junctions. At a frequency of 50 THz, we predicted tremendous loca- 
lisation of the electric-field intensity of the IDSPP to a modal area 
more than one million times smaller than the diffraction limit. As a 
new member of the existing plasmon family of bulk plasmons, sur- 
face plasmons, localised plasmons, etc., we foresee a new field of low- 
dimensional plasmonics based on IDSPPs, particularly branching 
out into various multi-material 2D systems^^"^^ 

Methods 

IDSPP Dispersion: Derivation and Analysis. We consider a hybrid 2D system in the 
y = 0 plane comprising two semi-infinite domains laterally connected along the z axis 
with sheet conductivities of a^^' = (7^^^' + ia^"^^" {x < 0) and a^^^ = a^^^' + /a^^^" {x > 0) 
(see inset in Fig. 1) immersed in a uniform medium with relative dielectric 
permittivity s. In the quasi-static theory of Volkov and Mikhailov^^'^^, based on the 
solution to Poisson's equation with an assumed electric potential of the form 
(p{r) = (po{x,y)exp{iq2Z—icQt), the dispersion relation of plasmons propagating along a 
ID junction between two adjoining 2D electron gases (2DEG) is given by^*^ 



1 , ^^^^ + u 
1+ ^ tanh 



-f- 



Here, Sapy = (T^^^py — o^^^py, where and y are tensor indices {x, z), and are the 
effective dielectric permittivities of the left (L) and right (R) 2DEG^^; 

= e + i^^xx^^ {co)q/(28QCo), where q = q^x + (fz is the permittivity of free 

space. For electrically isotropic 2DEGs in the absence of an external magnetic field, we 
have cr^^'^ xz = 0 and do^^ = 0. Taking doy^ = do = o"^^' — cr^^\ we obtain 



- ida coth 



71 J 1 



dc^ 



In 



£L(q,CL») 



When do 0, the dispersion relation reduces to when the hyperbolic cosine of the 
argument in the above equation is zero. To simplify the analysis we start by setting 
^(R,L)' _ -y^hich neglects ohmic and interband losses (the complex- conductivity case 
is treated in Supplementary Information), e.g., pristine graphene in the local and zero- 
temperature limit^° at a normalised frequency Q < 2; Q = /za»//i, where ji is the Fermi 
energy. Then, making the variable substitution sin0=(l + £,^)"^^^, we arrive at equa- 
tion (1) with the prescribed definitions. Inspection of equation (1) reveals that the left 
hand side of the plasmon dispersion can be written as the difference of two integrals 

.71/2 

j{KN) - f{N), where/(C) = d9 ln(sin 9-C). While the solution of this integral 
Jo 

lacks a simple closed-form expression for arbitrary C, its numerical solution com- 
bined with its asymptotic behaviour is revealing. We find that Im[/(C)] = 0 when 
C < 0, Im[/(C)] = 71 V2 when C > 1, and Im[/(C)] increases monotonically as C 
increases within the range 0 < C < 1 (see the inset of Fig. 4a; dashed red and blue 
curve). This behaviour immediately precludes solutions for purely real N (i.e., non- 
leaky) when i<C> 0 (i.e., metal-metal systems) because 0 <f{KN),f{N) < 71^/2, which 
constrains — 71^/2 < lm\f{KN) — f{N)] < so that the imaginary part of equation 
(1) could not be satisfied for any integer m. Indeed, this was earlier pointed out by 
Mikhailov^^ and corresponds to leakage of the so-called inter-edge plasmon into 
planar 2D plasmons towards the region with the smaller carrier density (i.e., smaller 
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Figure 4 | IDSPP dispersion and graphene conductivity, (a), The 
numerical solution [solid blue curve; equation (1)] and the asymptotic 
expression [dashed red curve; equation (2) for \ K\ > 0.5] are both shown. 
Inset: Real (solid curve) and imaginary (dashed curve) parts of the function 
/(C). The colour distinguishes C> 0 (red) from C < 0 (blue). The black 
overlaid dotted curve corresponds to Re[/(C + A)] ~ Re[/( — C)]. 
(b), Conductivity of graphene in the local and zero-temperature limit. 

Now, we turn to the pertinent case of the ID metal- dielectric (IDMD) interface [K 
< 0; a^^^" > 0, a^^^" < 0)] . This case has not been previously considered^*^'^^ because of 
the strictly positive sign of a" in conventional 2DEGs. New opportunities to realise 
IDMD interfaces in hybrid 2D systems prompted us to investigate fundamental 
solutions to equation (1) in 2D metal-dielectric systems. Our previous observations 
imply non-leaky solutions to equation ( 1 ) for K<0 when AT > 1 (and m = 0), provided 
that Re[/(- |K|iV)] = Re[/(iV)] (which satisfies the real part of equation (1)); indeed, 
lm\f{KN)] = 0 for iVK < 0 (i.e., iV > 0,K < 0), and Im[/(iV)] = nV2 for N > 1, thus 
satisf)^ing the imaginary parts of equation ( 1 ) . From the numerical integration of/(C), 
we empirically find the relation Re[/(C + A)] « Re[/(-C)] (C> 0) (see the inset of 
Fig. 4a; overlaid dotted black curve). Setting C = |K|iV leads to Re[/(|K|iV + A)] « 
Re[f(- |K|iV)], and recalling the solution condition Re[/(- |K|N)] = Re[/(iV)], we 
must have (comparing the positive arguments of the function Re[/(C)] which 
monotonically increases for C> 1, and noting A>l)|K|Ar+A = Ar. This indicates a 
dispersion relation of the form N = NJ{1 — \K\), where we have identified A = AT^ as 
the bare edge plasmon dispersion (i.e., for when K =0), or recalling the definitions of 
N and K and taking^*^ Ne^^/Vfl (the known approximate factor of the bare edge 
plasmon dispersion), we can write the following for the propagation constant: 

<j~2a/3/2££oCo/ [o^^^" + o^^^"^ . One may also consider the asymptotic behaviour in 

the limit 1; writing equation (1) as the difference of two integrals 

/( — \K\N) —f{N) and Taylor-expanding the respective arguments to the first order, we 

find 



1+T 



\K\JH1/\K\) 



; iV,N|K|> >1, 



(2) 



which further illustrates the limiting behaviour as \K\ ^ 1 of N ~ (2/7r)(l + |K|)/ 
(1 — |i<C|), which is consistent with the (1 — \K\)~^ dependence in the empirical 
relation obtained above. Solution to equation (1) only when m = 0 naturally pre- 
cludes the existence of higher-order, multipolar modes because of the step-like 
transition of the conductivity across the junction; this result is expected, considering 
that multipolar modes of a bare edge plasmon manifest only when the conductivity 
decreases monotonically from a constant to zero over a non-zero length a beyond 
some cutofP\ We note that the dispersion of the bare edge plasmon^*^'^^ is recovered in 

p/2 / ^ \ 

the limit that K = 0, given by the root of dO ln[ — — - — 1 ) = 0, which is satisfied 



by AT = ATg ~ 1.217. . In Fig. 4a, we plot the dispersion N{K) as determined from the 
numerical solution of equation (1) and as given by the analytical equation (2) given 
for |i<C| > 0.5 (see solid and dashed curves, respectively). 

Numerical analysis of IDSPP in hybrid graphene/graphene system. The IDSPP 
dispersion was obtained using the mode analysis tool in the COMSOL RF module 
(www.comsol.com). In our example of the graphene/graphene hybrid 2D system, the 
conductivities on either side of the junction are calculated using the random phase 
approximation (RPA) in the local and zero-temperature limit^° (Fig. 4b): 



0(Q(^'L)-2) 



-In 



2 + Q(^'L) 



2-Q* 



(3) 



In our COMSOL simulations, the graphene is incorporated into the numerical 
simulations as a thin film with a thickness of S = 0.2 nm, an effective relative 
dielectric permittivity^ of £^^'^^=l + iG^^'^V(jJSQS, and a^^'^^ given by equation (3). We 
note that the band structure of graphene is implicitly included in the permittivity 
through the conductivity as determined from the RPA. The almost perfect agreement 
between our numerical results and theoretical calculations (according to equation 
(1)) verif)^that the numerical results are sufficiently converged to the S^O limit (at S 
= 0.2 nm). 

From our simulations, we also observed the onset of plasmon leakage towards 
region R (x > 0) as (Q^^^^ >0.6, thus confirming the previous assertion^^ for the case 
when K > 0. At values of (Q^^^)~^ < 0.5, the dispersion must be solved for a complex 
value of ATbecause of the introduction of non-zero cr^^^'; nevertheless, considering that 
we would likely have a^^^' > (see Fig. 4b), the IDSPP will be strongly damped. 
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